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Abstract 

We derive improved upper bounds on the number of crossing-free straight-edge spanning 
cycles (also known as Hamiltonian tours and simple polygonizations) that can be embedded 
over any specific set of N points in the plane. More specifically, we bound the ratio between the 
number of spanning cycles (or perfect matchings) that can be embedded over a point set and the 
number of triangulations that can be embedded over it. The respective bounds are 0(1.8181''^) 
for cycles and 0(1.1067^) for matchings. These imply a new upper bound of 0(54.543^) on the 
number of crossing-free straight-edge spanning cycles that can be embedded over any specific 
set of N points in the plane (improving upon the previous best upper bound 0(68.664^)). Our 
analysis is based on Kasteleyn's linear algebra technique. 

1 Introduction 

Finding a Hamiltonian cycle in a graph is NP-complete even if the graph is known to be planar 
[7]. Counting the number of Hamiltonian cycles that are contained in a graph is #P-complete 
even if the graph is known to be planar and all of its vertices have degree three |15] . In this paper 
we consider a more general case (and thus, probably more difficult) concerning the number of all 
crossing-free straight-edge Hamiltonian cycles that can be embedded over a specific set of points 
in the plane. That is, given a set S of labeled points in the plane, we consider the number of 
Hamiltonian cycles that have a straight-edge planar embedding over S. 

We only consider the problem of obtaining an upper bound for the number crossing-free straight- 
edge Hamiltonian cycles that can be embedded over a set of N points in the plane. To do this, 
we rely on Kasteleyn's technique [H], and on edge-flipping techniques that were developed in a 
previous paper by the authors [9]. No familiarity with [9] is necessary, since we re-introduce all the 
notions that we require from it. We now give a detailed and more formal definition of the problem. 

A planar graph is a graph that can be embedded in the plane in such a way that its vertices 
are embedded as points and its edges are embedded as Jordan arcs that connect the respective 
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pairs of points and can meet only at a common endpoint. A crossing-free straight-edge graph is a 
plane embedding of a planar graph such that its edges are embedded as non-crossing straight line 
segments. In this paper, we only consider crossing-free straight-edge graphs. Moreover, we only 
consider embeddings where the points are in general position, that is, where no three points are 
collinear. (For upper bounds on the number of graphs, this involves no loss of generality, because 
the number of graphs can only grow when a degenerate point set is slightly perturbed into general 
position.) For simplicity, we sometimes refer to such graphs as plane graphs. 

We focus on upper bounding the maximal number of plane spanning cycles (also known as 
Hamiltonian cycles, Hamiltonian tours, and simple polygonizations) that can be embedded over 
a fixed set of points in the plane. For a set S of points in the plane, we denote by C{S) the set 
of all crossing-free straight-edge spanning cycles of 5, and put sc(5) := |C(5)|. Moreover, we let 
5c{N) = max|5|^^ 50(5). So, in other words, the main goal of this paper is to obtain sharp upper 
bounds on sc(A^). 

There are many similar variants of this problem, such as bounding the number of plane forests, 
spanning trees, triangulations, and general plane graphs. Recent work on some of these variants can 
be found in [HO [21], and we try to keep a comprehensive list of the up-to-date bounds in a dedicated 
webpag^. It seems that the case of spanning cycles is the most popular one, already considered in 
[21 [21 m O [6l [T71 [21] and many others. Moreover, spanning cycles were the first case for which bounds 
were published, namely the bounds 3/20 -10^/^ < sc(A^) < 2-6^-^ ■ {[N/2\)\ in [17]. A brief history 
of the steady progress on bounding the number of spanning cycles can be found in a dedicated 
webpage by Erik Demain^. Currently, the best known lower bound is sc{N) = 0(4.642^), due to 
Garcia, Noy, and Tejel [6], and the previous upper bound is sc(A^) = 0(68.664''^) by Dumitrescu et 
al. [5]. We derive the improved bound sc(A^) = 0(54.543^). 

These problems have also been studied from an algorithmic point of view, deriving algorithms 
for enumeration or counting of the plane graphs (or other graph types) that can be embedded over 
a given point set (such as in [131 119j). The combinatorial upper bounds are useful for analyzing 
the running times of such algorithms, and also to answer questions such as "how many bits are 
required to represent a triangulation (or any other kind of plane graphs)?". 

Our bound (as do some of the previous bounds) relies on triangulations. A triangulation of a 
set S oi N points in the plane is a maximal plane graph on S (that is, no additional straight edges 
can be inserted without crossing some of the existing edges). For a set S of points in the plane, 
we denote by T{S) the set of all triangulations of S, and put tr(5) := |T(S')|. Moreover, we let 
tr(A^) = max|5|=7v ti'('S')- Currently, the best known bounds for tr(A^) are tr(A^) < 30^ [21], and 
tr(iV) = n{8.65^) [5]. 

The upper bound by Dumitrescu et al. [5] is obtained by proving that for every set S of 
points in the plane 5c{S) = O (2.2888^) • tr(S'). This has sharpened an earlier bound of Buchin et 
al. [1|, who showed that every triangulation T of 5 contains at most 30^/^ ~ 2.3404^ spanning 
cycles (i.e., cycles whose edges belong to T), implying that sc{S) < 2.3404^ • tr{S). Combining 
the above ratio with the bound tr(A^) < 30^ directly implies the asserted bound. We derive our 
bound in a similar manner, showing that sc(5) = O (1.8181^) • tr(S) = 0(54.5430^). 

In spite of our improved bound, we strongly believe, and conjecture, that for every point set 



http: //www. cs . tau. ac . il/~ shef f era/ counting/PlaneGraphs .html [(version of June 2010). 



http://erikdemaine.org/polygonization/ (version of June 2010) 



The implication comes from the fact that every spanning cycle, and in fact every plane graph, is contained in at 
least one triangulation; see Section (2] 
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Figure 1: Two spanning cycles embedded over a double chain point configuration. 



S (of size at least Nq, for some constant A'^o) one has sc(S') < tr(S'), and perhaps even a much 
sharper ratio holds. The best lower bound for this ratio that we know of is obtained from the 
double chain configuration, presented in (Uj (and depicted in Figure [T]) . It is shown in [B] that when 
S" is a double chain configuration, tr(5) = Q* (8^) and sc(5) = (4.64^)0 Thus, in this case, 
sc{S)/tr{S) = n*{0.58'^). (It is stated in [I], albeit without proof, that sc(5) = 0(5.61^), so this 
example supports our conjecture.) 

In Section [2] we go over the preliminaries required for our analysis. These include, among 
others, the edge-flip techniques used in [9]. Section [3] derives the bound 50(5") = O (l2^/^) ■tr(5) = 
O (l.8613^) • tr(S') for any set S of N points in the plane. As part of this derivation, we describe 
Kasteleyn's technique for counting perfect matchings and present a new way of applying it. In 
Section |4] we use the same methods to prove an upper bound on the ratio between the number of 
plane perfect matchings and the number of triangulations, showing that pm(S') = 0(1.1067^) •tr(S') 
(where pm(S') is the number of matchings that can be embedded over the point set S). Finally, 
Section [5] contains a more complex analysis of spanning cycles, implying the improved bound 
sc{S) = O (10.9247^/4) • tr(5) = O (1.8181^) • tr{S). 

2 Preliminaries 

In this section we establish some notations and lemmas that are required for the following sections. 

Given two plane graphs G and H over the same point set S, if every edge of G is also an edge of 
H, we write G H. 

Hull edges and vertices (resp., interior edges and vertices) of a graph embedded on a point set S 
are those that are part of the boundary of the convex hull of S (resp., not part of the convex hull 
boundary). 

Given a set S of N points in the plane, we denote by h the number of hull vertices of S, and put 
n = N — h, which is the number of vertices in S interior to its convex hull. 

2.1 The support of a graph 

Let us denote by sc^(A^) the maximal number of plane spanning cycles that can be contained in 
any fixed triangulation of a set of N points in the plane. Moreover, denote the set of spanning 
cycles contained in a triangulation T by C{T), so 5C^{N) = max|5|=7v, T€T{S) 1^(^)1- 

Any spanning cycle (or, for that matter, any plane graph) is contained in at least one triangu- 
lation. Therefore, we can upper bound the number of spanning cycles of a set S of N points in 
the plane by going over every triangulation T G T{S) and counting the number of spanning cycles 

*In the notations 0*(), 0*(), and we neglect polynomial factors. 
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contained in T. This implies the bound sc(A^) < tr(A^) •sc^(A^). Applying the bounds tr(A^) < 30^ 
from [21] and sc^(A^) < 30^/'' from [1], we obtain sc{N) < 30^^/"^ ^ 70.21^. 

This bounding method seems rather weak since it potentially counts some spanning cycles many 
times. For example, consider a spanning cycle consisting of two convex chains facing each other, 
as depicted in the left-hand side of Figured] Garcia, Noy, and Tejel [6] show that such a spanning 
cycle is contained in 0*(8^) triangulations of its point set. Therefore, the above method will count 
this spanning cycle 0*(8^) times. However, as stated in [1], this point set has only 0(5.61^) 
spanning cycles. 

In order to deal with this inefficiency, we define the notion of support (the same notion was 
also used in [U El EH [22l [25]). Given a plane edge graph G embedded over a set S of points in the 
plane, we say that G has a support of x if G is contained in (exactly) x triangulations of S; we 
write supp(G') = x. Notice that 



because every spanning cycle G contributes exactly one to the right side of the equation (it appears 
in supp(C) terms of the first sum, and contributes l/supp(C) in every appearance). We will use 
([1]) to obtain better upper bounds for sc(iV), by showing that, on average, supp(C) is large. 

2.2 Ps-flippable edges 



Figure 2: (a) The edge ce can be flipped to the edge ad. (b) A set of three (dashed) ps-flippable edges which 
are diagonals of interior-disjoint convex quadrilateral and convex pentagon, (c) A convex decomposition obtained 
by removing the dashed edges from the triangulation; here too the dashed edges are diagonals of pairwise disjoint 
convex faces, and form a set of ps-flippable edges. 

An edge in a triangulation is said to be flippable, if its two incident triangles form a convex 
quadrilateral. A flippable edge can be flipped, that is, removed from the graph of the triangulation 
and replaced by the other diagonal of the corresponding quadrilateral. Such an operation is depicted 
in Figure [2|^a) , where the edge ce can be flipped to the edge ad. 

In [9], we present the concept of pseudo simultaneously flippable edges (or ps-flippable edges, for 
short). Given a triangulation T, we say that a subset F of its edges is a set of ps-flippable edges if 
the edges of F are diagonals of interior-disjoint convex polygons (whose boundaries are also parts 
of T). For example, in Figure [5]^b) , the three dashed edges form a set of ps-flippable edges, since 
they are diagonals of interior-disjoint convex quadrilateral and convex pentagon (another set of 
ps-flippable edges, in a different triangulation, is depicted in Figure El^c)). 

Ps-flippable edges are related to convex decompositions. A convex decomposition of a point set 
5 is a crossing- free straight-edge graph D on S such that (i) D includes all the hull edges, (ii) each 
bounded face of is a convex polygon, and (iii) no point of S is isolated in D. See Figure [2|^c) 




(1) 




(a) 



(b) 
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for an illustration. For additional information about convex decompositions, see, for example, |10j . 
Notice that if T is a triangulation that contains -D, the edges T \ D form a set of ps-flippable 
edges, since they are the diagonals of the interior-disjoint convex polygons of D (again, consider 
the dashed edges in Figure [2^c) for an illustration). Thus, finding a large set of ps-flippable edges 
in a triangulation T is equivalent to finding a convex decomposition with a small number of faces 
(or edges) in T. 

In [9], we prove the two following lemmas. 

Lemma 2.1 Every triangulation T over a set of N points in the plane contains a set F of N/2 — 2 
ps-fiippbale edges. Also, there are triangulations with no larger sets of ps-fiippable edges. 

Lemma 2.2 Consider a triangulation T, a set F of N/2 — 2 ps-flippable edges in T, and a graph 
G (^T. If G does not contain j edges from F then supp(G) > 2^ . 

Proof sketch. Consider the set F' = F \ G oi j ps-flippable edges. The convex faces of T \ F' 
can be triangulated in at least 2^ ways (the bound is tight when every edge of F' is a diagonal of a 
distinct quadrangular face of T \ F'), and each of the resulting triangulations contains G. See [9] 
for more details. 

□ 

(b) 

Figure 3: (a) A vertex is valid if and only if it is not a reflex vertex of any face, (b) In a point set of an even 
size, every spanning cycle is the union of two edge-disjoint perfect matchings. 

We now describe another property of convex decompositions (not discussed in [9]). Consider 
a set S of points in the plane and a crossing-free straight-edge graph G embedded on S. We say 
that an interior point p £ S has a valid triple of edges in G if there exist three points a,b,c G S 
such that p is contained in the convex hull of {a, b, c} and the edges ap, bp, and cp belong to G. 
To simplify the notation, we refer to vertices with valid triples as valid (with respect to G), and to 
the other interior vertices as non-valid. 

Lemma 2.3 Let S be a set of points in the plane and let G be a crossing-free straight-edge graph 
over S that contains all the edges of the convex hull of S. Then G is a convex decomposition of S 
if and only if every interior vertex of S is valid with respect to G. 

Proof. An interior vertex v is a reflex vertex of some face of G if and only if v is non-valid 
(an example is depicted in Figure EJ^a)). The lemma follows by observing that G is a convex 
decomposition if and only if no bounded face of G has a reflex vertex. □ 

2.3 Spanning cycles and perfect matchings 

Our analysis, as most of the previous works dealing with the number spanning cycles, heavily relies 
on the number of plane perfect matchings on S (for example, see [ll[5l[23]). To see the connection 




(a) 
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between the two problems, notice that if 15"! is even, every spanning cycle C is the union of two 
edge-disjoint perfect matchings on S; namely, the matching consisting of the even-indexed edges 
of C, and the matching consisting of the odd-indexed edges. An illustration of this property is 
depicted in Figure Ell^b). Denote by M{S) the set of all plane perfect matchings on 5, and put 
pm(5') = |A^(S')|. We also set pm(A^) = max|5|=jv prn('S')- Hence, a simple upper bound on sc(5) 
is pm(S')^. In general, the union of two edge-disjoint perfect matchings is not always a spanning 
cycle, but it is a cover of S by vertex-disjoint even-sized cycles. 
To deal with point sets of odd size, we use the following lemma: 

Lemma 2.4 Let c > 1 be a constant such that every set S of an even number of points in the plane 
satisfies sc(5) = 0{c^^^). Then sc{S) = 0{S^^) also holds for sets S of an odd number of points. 

Proof. Consider a set S oi N points in the plane, where N is odd. Pick a new point p outside the 
convex hull of S, and put S' = S U {p}. Let C be a plane spanning cycle of S. Then there exists 
an edge e = vu of C such that p can be connected to the two endpoints u,v oi e without crossing 
C. Indeed, this is a projective variant of the property, noted in [8], that every finite collection of 
non-crossing straight segments in the plane contains a segment e such that no other segment lies 
vertically above any point of e (see also \18\ Section 8.7]). By replacing e with the edges vp and 
pu, we obtain a crossing-free spanning cycle of S' . This implies that we can map every spanning 
cycle of S" to a distinct spanning cycle of S' , and thus, sc{S) < sc{S'). The lemma then follows 
since sc(5') = 0{c^+^) = O(c^). H 

Bounding the number of perfect matchings on S within a fixed triangulation T can be done by 
the beautiful linear-algebra technique of Kasteleyn [13], described in detail in [161 Section 8.3]; see 
Section [3] for more details. Buchin et al. [3] have used this technique to show that any triangulation 
T of 5" contains at most 6^/^ perfect matchings, and at most 30^/'^ ~ 2.3403^ spanning cycles. 
We also note that Sharir and Welzl [24] showed that pm(S') = O (l0.05^), completely bypassing 
the approach of counting matchings (or other graphs) within a triangulation. 



3 A first bound 

In this section we first review an enhanced variant of Kasteleyn's technique and then use it to derive 
the following initial upper bound on sc{S). 

Theorem 3.1 For any set S of N points in the plane, 

sc{S) = O (l2^/^) • tr(5) = O (l.8613^) • tr{S). 



Proof. First, by Lemma l2.4( we may assume that is even. Consider a triangulation T of S. 
As already observed, every spanning cycle contained in T is the union of two edge-disjoint perfect 
matchings contained in T. Given a plane graph G, we denote by Ai{G) the set of all perfect 
matchings that are contained in G. Recalling ([T]), we have 



supp(Mi U M2) 

AI^,Al2 edge-disjoint 

(The inequality comes from the fact that not every pair Mi,M2 of matchings, as in the sum, 
necessarily yields a spanning cycle.) Let us fix the "first" perfect matching Mi C T; as mentioned 
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(a) 



(b) 



Figure 4: (a) The set {e, ei, 62} is a valid triple of edges, (b) A set of 12 points in a double circle configuration. 

above, Buchin et al. [1] prove that |A^(r)| < 6^/^, so there are at most 6^/^ choices of Mi. Next, 
we construct a convex decomposition D such that Mi C D C T, as follows. We start with Mi and 
add all the missing hull edges; let us denote the resulting graph as D' . By Lemma |2.3| it suffices 
to ensure that every interior point p & S is connected in D to (at least) three points a,b,c € S, 
such that p is inside the convex hull of {a, b, c}. Every interior vertex p oi S has degree 1 in D' , so 
we start by setting D := D', and then, for each interior point p € S", we add to D two additional 
edges of T adjacent to p, so as to create a valid triple. To do so, let e be the edge of D' (that is, of 
Ml) incident to p, and let A be the ray emanating from p in the opposite direction. Let ei (resp., 
62) be the first edge of T incident to p encountered in clockwise (resp., counterclockwise) direction 
from A; see Figure ID^a). Then {6,61,62} is a valid triple of edges, and we add 61,62 to D. After 
applying this step to each interior point p, the resulting graph D is indeed a convex decomposition 



We denote by F the set of edges that are in T but not in D. The edges of F are diagonals 
of interior-disjoint convex polygons, and thus F is a set of ps-flippable edges. By Euler's formula, 
the triangulation T contains 3A^ — 2h — 3 interior edges, and D contains at most 2n + N/2 interior 
edges (at most N/2 edges of Mi and 2n added edges to form n valid triples). Therefore, 



Remark. Note the strength of this bound: Lemma l2.ll has a rather involved proof, given in j9], 
and it yields a set of N/2 — 2 ps-flippable edges in the entire triangulation. In contrast, here we 
get the same number (minus 1) after we remove from T an arbitrary perfect matching, with a 
considerably simpler analysis. Thus the significance of the analysis in [9] (involving Lemma 12. ip 
is only for triangulations which contain no perfect matching on S. For example, any triangulation 
with more than N/2 interior vertices of degree 3 cannot contain a perfect matching, since, as is 
easily checked, two interior vertices of degree 3 cannot share an edge. 

Without loss of generality, we assume that F consists of exactly N/2 — 3 edges. We now proceed 
to bound the number of ways to choose the second matching M2 while taking the supports of the 
resulting graphs Mi U M2 into account. Since Mi and M2 have to be edge-disjoint, we can remove 
the N/2 edges of Mi from T, and remain with a subgraph T' that has fewer than 5N/2 edges. 
Next, we define a weight function /_f over the edges of T', such that every edge in F has a weight 
of 1 and every other edge has a weight of 1/2. We define the weight ^(M2) of a perfect matching 
M2 C T' as the product of the weights of its edges. Therefore, if M2 contains exactly j edges of 
F, then At(M2) = (l/2)^/2-j\ Moreover, for such a matching M2, we have \F \ M2I = N/2 -3-j. 



of S. 



F\>3N-2h-3-{2n + N/2) = N/2 - 3. 
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Clearly, F \ M2 is also a set of ps-flippable edges, none of which belongs to Mi U M2. We thus have 

1 1 



< 



8^(Af2 



supp(Mi U M2) - 2^/2-3-i 
which implies that, given a specific triangulation T and a specific perfect matching Mi C T, 

1 



E 



with T' = T\ Ml, as above. 



supp(Mi U M2) 



< 8 ^(^^2) 



(2) 



Af2e>!(T') 



Kasteleyn's technique: An enhanced version. We now apply an extension of Kasteleyn's 
technique to estimate the sum in the right-hand side of ([2]) . Here is a brief overview of the technique 
being used (where instead of the original technique, we apply a weighted extension of it). Given 
an oriented graplil G = {V, E) with no anti-parallel edges and a weight function ji over the edges, 
we define the following weighted adjacency matrix ^ = {bij)iyxN of {G,ii), 

fi{e), if e = G E, 
-/i(e), if e = G E, 
0, otherwise 

(where N = \V\, and the rows and columns of ^ correspond to an arbitrary fixed enumeration 
of the vertices) . 

An easy extension of Kasteleyn's theorem states that every planar graph G can be oriented into 
some digraph G such that, for any real- valued weight function fi on its edges, we have 

2 



E ^(^^) = det J (3) 

(recall that /u(M) = neGA//"(^))- "pure" form of Kasteleyn's theorem = 1 (i.e., G 

is unweighted) and the left-hand side is just the squared number of perfect matchings in G. A 
detailed presentation of Kasteleyn's theorem can be found in [161 Section 8.3]. The extension to 
weighted graphs is given in Exercise 8.3.9 therein. 

We denote by hi the column vectors of B, for 1 < i < N , and estimate the above determinant 
using Hadamard's inequality 



det B, 



N 

<m"' 

i=l 



'i 2- 



(4) 



Applying the above machinery to our plane graph T' (i.e., using ([3]) and (j4])), with the edge weights 
fi as defined above, we have 

\ 1/4 



M2eM(T') 



det(5^, 



< 



N 



N 



\i=l 



< 



1 ^ 

-TWh 



N/4 



N 



N/4 



(5) 



egT' 



'''We follow here the notation used in [16] to denote a digraph obtained from an underlying undirected graph by 
giving each of its edges an orientation. 
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(where we have used the arithmetic-geometric mean inequahty and the fact that every edge of T' 
has two corresponding matrix entries). We note that the bound Q^/^ on the number of perfect 
matchings in a triangulation T is obtained in [5] by applying the unweighted version of Kasteleyn's 
theorem to the entire T. In this case ^eeT /^(^)^ is the number of edges of T, which is at most 3A'', 
and the bound fohows. 
By noting that 

\T' \F\< 5N/2 - {N/2 - 3) = 2iV + 3, 
and combining this with ([2]) and ([5]), we obtain 

N/4 



M2£MiT') ' \ eeT' / 



< 8- -• l^-(7V/2-3)+(-) •(2Ar + 3) 



O (2^/4) . (6) 



Recalhng once again that a triangulation contains at most 6^/^ perfect matchings [3] (that is, there 
are 6^/^ ways of choosing Mi), and combining this with ([6|), we obtain 

sc(S)< y y -^-^< y 6^/4.0(2^/4) =ofl2^/4) .tr(5), 

^ ^ - ^ ^ supp(AfiUM2 ~ ^ \ ) V ; ^ ^' 

T&T{S) Mi,M2eM(T) f^^^ ^ Ter(5) 

Af]^,A/2 edge-disjoint 

as asserted. H 

Combining Theorem 13.11 with the bound tr(A^) < 30"^ [21], we obtain 
Corollary 3.2 sc(iV) = O (55.8363^). 



4 Perfect matchings and triangulations 

In this section we apply the machinery of the previous section to derive an upper bound on the 
ratio between the number of plane perfect matchings and the number of triangulations. As already 
mentioned in Section [2l Kasteleyn's technique implies that a triangulation of a set of points can 
contain at most 6^/4 perfect matchings (see [1]). This implies that for every set S of N points in 
the plane, pm(5) < 6^/4 • tr(5) « 1.5651^ • tr{S). We will improve this bound, using lower bounds 
on the supports of perfect matchings, in a manner similar to that in Section [3j 

Before proceeding, we note the following lower bound on the ratio pm(5')/tr(S'). Let S* be a 
double circle configuration, depicted in Figure [D^b), consisting of N points (see [11] for a precise 
definition). An inclusion-exclusion argument implies that tr(S') = 12^/^ (see |1H I20j). Moreover, 
Aichholzer et al. [Tj proved that pm(S') = ©*(2.2'^). Therefore, in this case, pm(S')/tr(S') f« 
e*(0.635^). 

We now present an improved upper bound for this ratio. 
Theorem 4.1 For any set S of N points in the plane, 

pm(5) < 8 • (3/2)^/4 . tr{S) = 0(1.1067^) • tr(5). 
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Proof. The exact value of pm(S') is 



pm(5)= y y — . (7) 

^ ^ supp(M) ^ ' 

Consider a triangulation T € 7'(S') and a perfect matching M ^ T. As shown in the proof of 
Theorem 13.11 there exists a set of N/2 — 3 ps-flippable edges in T \ M. Therefore, the support of 
M is at least 2^/^~^. Combining this with ([7]) implies 

pm(5)< E E :^ = 8-(3/2)^/'-M5)- 



Ter(S) MeM{T) TeT{S) 



As already mentioned above, this does not imply a new bound on pm(A^), since Sharir and 
Welzl [23] showed that pm(S') = O (l0.05^), bypassing the approach of counting matchings within 
a triangulation. We are not aware of any construction for which pm(S') > tr(S'), and offer the 
conjecture that there exists a constant c < 1 such that pm(5) = 0{c^^^ ■ tr{S)) for every finite 
set S of points in the plane. (See also the conjecture concerning spanning cycles, made in the 
introduction.) 



5 An improved bound 

In this section we present a more complex analysis for the number of spanning cycles, obtaining 
a slightly better bound than the one presented in Section [3l The analysis has three parts, each 
presented in a separate subsection. 

Let us denote the number of interior vertices of degree 3 in the triangulation T as f3(T'). 
Moreover, let us denote the number of flippable edges in T as flip(T). In Subsection 15.11 we give 
an upper bound for J2ceC{T) supp(C) ^^^^ depends on f3(T). In Subsection 15.21 we give an upper 
bound for '^c'eC{T) supp(c) depends on flip(T). Finally, in Subsection 15.31 we combine these 
two bounds to obtain 

sc(iV) = 0(10.9247^/'') • tr(iV) = 0(1.8181^) • tr{N) = 0(54.5430^). 



5.1 A t;3(T)-sensitive bound 

In this subsection we derive the following bound, which is a function of N and U3(T). 

Lemma 5.1 Let T be a triangulation over a set S of N > 6 points in the plane, such that N is 
even and S has a triangular convex hull; also, let V3{T) = tN . Then 

V 1 ^J^ ( {2-t){2-t/2) y^'Y^' 

Proof. As before, we treat every spanning cycle as the union of two edge-disjoint perfect matchings 
Mi,M2 G Ai{T). We start by bounding the number of ways to choose the first perfect matching 
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Ml. For this, we use the standard variant of Kasteleyn's technique, with the weight function fj, = 1 
(i.e., G is unweighted). 

RecaU the inequahty X^jv/g_A4(T) A^-^) ^ (lli^i ll^illi) obtained in Equation ([5]), where the 

bi (for 1 < i < N) are the column vectors of the adjacency matrix of the oriented graph T. 
Substituting /i = 1, the left hand side becomes the number of perfect matchings in T, and the 
squared norm of each column vector is the degree of the vertex corresponding to that column. 
Since every column that corresponds to a vertex of degree 3 has a squared norm of 3, the product 
of the squared norms of these is 3"3{t) ^ ^tN ^ 

For the remaining N — V3{T) columns, we use, as in Section [3l the arithmetic-geometric mean 
inequality to bound the product of their squared norms (as in Equation ([5])). This yields the bound 

(7V_„3(T))/4 / ^ N(JV{l-t))/4 

(8) 



N-V3{T)J \N{l-t) 

where X is the sum of the degrees of all vertices other than those counted in vsiT). The sum of 
the degrees over the vertices of any specific triangulation is smaller than 6A^, and the sum of the 
degrees of the interior degree-3 vertices in T is 3v3{T). Therefore, we have 

X <QN -3vs{T) = 3N{2-t). (9) 

Combining ([8]), ([9]), and the product of the squared norms that correspond to interior vertices of 
degree 3, implies that the number of ways to choose Mi is less than 

(10) 

Next, let us fix a specific perfect matching Mi G Ai{T). As shown in the beginning of the proof of 
Theorem 13. H there exists a set F of N/2 — 3 ps-flippable edges in T, none of which belongs to Mi. 

We continue as in the proof of Theorem 13. H by assigning a weight of 1 to the edges of F and 
a weight of 1/2 to the rest of the edges of T \ Mi, and then applying Kaseteleyn's technique to 
bound the sum 




N 



y — — r < 8 y u(M2) < 8 ( TT \\b'\ 

M2£M{T') AheMiT') \i=l 



1/4 



where 5^ are the column vectors of the oriented adjacency matrix of T \ Mi (recall ([2])). 

An interior vertex v of degree 3 in T has only two edges adjacent to it in T \ Mi, both not in 
F (since an edge adjacent to an interior vertex of degree 3 cannot be flippable). Therefore, the 
squared norm of a matrix column that corresponds to such a vertex is (1/2)^ + (1/2)^ = 1/2, and 
the product of the squared norms of all such columns is 1/2^^^^^) = 1/2*^. 

For the remaining — V3{T) columns, we may once again use the arithmetic-geometric mean 
inequality to obtain a bound similar to the one in ([8]). However, this time we get a different value 
for X, since (i) some of the edges of T were removed, and (ii) some of the remaining edges were 
reweighted. The edges of F have remained and still have a weight of 1 each, so they contribute at 
most 2 • (A^/2 — 3) • 1 < to X. Every other edge contributes 2 • 1/4 = 1/2 if it is not incident to 
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an interior vertex of degree 3 in T, and 1/4 otherwise. Since a triangulation has fewer than 3A^ — 3 
edges, there are fewer than 2N edges in T \ {F U Mi}, and we get 

X<N + {2N- 2v3{T)) ■ ^ + 2v3{T) ■ ^ = 2N - = iV(2 - t/2). 

By combining this with the rest of the squared norms and with the present version of ([8]) , we have 

Finally, to complete the proof, we combine (fTOl) and (fTTI) . and obtain 



< 3 



y y I 

M-^,M2 edge-disjoint 

^ Ar/4 / ^ ^ Ar/4 / ^ s i A 

2-A \ n /2-i/2V \ / 3 /(2-t)(2-t/2)^ 



1-ty / \2t \ 1 - 1 J / l2*V (1-0 



□ 



Remark. Notice that in the worst case (i.e., when t = 0) we obtain the same asymptotic value 
as in our initial bound of 12^/^. Similarly, the bound in (jlOp becomes 6^/'' when t = 0, as m 
Buchin et al. [4j. 

5.2 A flip(T)-sensitive bound 

Hurtado, Noy, and Urrutia [12] prove that that flip(T) > N/2 — 2, and that this bound is tight in 
the worst case (the upper bound is also implied by Lemma l2.lt see also [9]). In this subsection we 
obtain a bound as a function of flip(T), which improves our initial bound of 12^/^ when flip(T) is 
larger than N/2 by some positive fraction of N. We start by mentioning two basic properties of 
triangulated polygons. 



(a) (b) 



Figure 5: (a) The shaded faces are the two ears of the triangulated polygon, (b) A polygon with three flippable 
diagonals, two of which form a ps-flippable set, and with 5 = 2^(5/4)^ triangulations. 



Ears. Given a triangulated polygon P, an ear of P is a bounded face of the triangulation with 
two of its edges on the boundary of P. It can easily be shown that every triangulated simple 
polygon with at least four edges contains at least two ears (whose boundary edges are all distinct). 
For example, the two ears of the triangulated polygon in Figure [5]^a) are shown shaded. 
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1 f2N\ 

Catalan numbers. The A^-th Catalan number is Cm = . It is well known that a 

N + 1\N J 

convex polygon with N vertices has Cm-2 triangulations (e.g., see [23^ Section 5.3]). Therefore, 
the number of triangulations of a convex polygon with d > 1 diagonals is 

1 f2d + 2\ , ff>Y~^ , , 

(the inequality can easily be verified by induction). 

Next, we define Cgon as the maximum real number satisfying the following property. Every 
simple polygon P that has a triangulation Tp with k of its diagonals flippable and with I < k 
of these diagonals forming a ps-flippable set, has at least 2'cg(^' triangulations. Notice that the 
triangulations under consideration, including Tp, are triangulations of the polygon P, and not of 
its vertex set. 

5 

Lemma 5.2 x < Cgon ^ ^ with x ~ 1.17965 the unique real root of the polynomial 1 + 4x^ — Ax^. 

That is, every simple polygon P that has a triangulation Tp with k of its diagonals flippable and 
with I < k of these diagonals forming a ps-flippable set, has at least 2^x^~^ triangulations. 

Proof. Figure [5Kb) depicts a polygon that implies the upper bound. 

We prove the lower bound by induction on I and k. To have some base case for this induction, 
notice that when I = k, P has at least 2'^ = 2^x^~^ triangulations. 

Next, consider a polygon P and a triangulation Tp of P, such that Tp contains a maximal 
set F of I ps-flippable edges, and k — I > Q flippable edges not in F. Let /i,/2,-'' -.fj be the 
non-triangular faces of the convex subdivision Tp \ F. (Notice that Tp \ F is a convex subdivision 
of a polygon, and not of the convex hull of a point set.) Each /j is a convex polygon with mj > 4 
sides and — 3 diagonals, and thus, I = X^^=i(?Tij — 3). 

A flippable edge not in F must be on the boundary of some /j, since otherwise F is not maximal. 
We say that such a flippable edge e ^ F is covered by a face fi if (i) e is on the boundary of /j, 
and (ii) the triangle A in Tp that is contained in fi and incident to e has two edges in F (so e is 
the only edge of A on the boundary of fi). Ears are incident to two edges of their containing fi, 
and thus, cannot cover any edge. Since any triangulated polygon contains at least two ears, with 
four distinct boundary edges, a polygon with m sides (and m — 3 diagonals) can cover at most 
m — 4 edges. Therefore, if all flippable edges not in F are covered, then k — I < '^l^iimi — 4). By 
multiplying the number of triangulations of the /j's and applying ()12p . we get that the number of 
triangulations of P that contain Tp \ F is at least 

We are left with the case where there is a flippable edge e ^ F that is not covered by any fi . Let 
e' be the edge obtained by flipping e. We now derive a lower bound on the number of triangulations 
of P that contain e and on the number of triangulations of P that contain e'. To bound the number 
of triangulations that contain e, we partition P into two interior disjoint simple polygons P', P" 
by "cutting" P at e. More precisely, we consider the two triangulated polygons Tpi,Tpii C Tp. 
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(a) (b) 

Figure 6: (a) Removing the two (shaded) triangles incident to e' produces at most four triangulated sub-polygons, 
(b) A spanning cycle partitions the convex hull of the point set into interior-disjoint polygons. 

Together, these two triangulated polygons contain k — 1 diagonals that are flippable. Moreover, 
the set F remains a set of I ps-flippable edges which are diagonals of the two polygons. Thus, the 
induction hypothesis implies that there are at least 2'a;'^~^~' triangulations of P that contain e. 
To obtain a similar bound for the number of triangulations containing e', we produce at most four 
triangulated sub-polygons of Tp by removing the two triangles incident to e (which cover the same 
quadrilateral covered by the two triangles incident to e' in the new triangulation) . Such a case is 
illustrated in Figure [6l^a), where the triangles incident to e' are shaded. This partitioning may cancel 
the flippability of at most five edges (those incident to the two triangles adjacent to e). At most two 
out of the five edges may be ps-flippable, since e can only be incident to ears of polygons of Tp \ F. 
Using the induction hypothesis again, we get that there are at least 2^~'^ x^''~^'^~^^~'^'^ = 2^~'^x^~^~^ 
triangulations that contain e' (and the four edges around it). Therefore, P has at least 



2'x"~"~' + 2'~"x"~'"'^ = Tx"-' - + — = I'x"-' ^ = 2'x 

43.3 J 43.3 

triangulations, as asserted (recall that x is the root of the polynomial 1 -|- — 4x^). H 
Next, we show how to use Cgon and flip(T) to bound J2ceC{T) supp(C) • 

Lemma 5.3 Consider a triangulation T with flip(T) = N/2 — 3 -|- kN , for some k > 0, and let x 
be the constant presented in Lemma \5.S\ Then 

J2 < . A3 + (7' - + l/2))(4 + (x2 - ^ 

c 

where 



supp(C) ^ " V x"^" 



= X ■ e 4(4+(3:-!-1)k) _ 

Proof. Once again, we treat every spanning cycle as the union of a pair of edge-disjoint perfect 
matchings Mi,M2 G 7W(T), and use Kasteleyn's technique (as presented in Section [3]) to bound 
the number of such pairs. We start by fixing some perfect matching Mi € M{T) and denote the 
number of flippable edges of T that are in Mi as flipy(Mi). As shown in the proof of Theorem 13. H 
there is a set of at least N/ 2 — 3 ps-flippable edges in T \ Mi . We restrict our attention to a set F 
of exactly N/2 — 3 ps-flippable edges in T \ Mi. 

For the choice of M2, we define a weight function /i(-) on the edges of T \ Mi, such that 

2, if e e F, 
fj,{e) = ^ X, if e ^ F is flippable, 
1, if e is not flippable. 
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Notice that any spanning cycle partitions the convex hull of its point set into interior-disjoint simple 
polygons; an example of such a partition is illustrated in Figure [5{b). The support of the spanning 
cycle is the product of the number of triangulations of each of these polygons. For a fixed choice 
of M2 (and of Mi), denote by Pi, ... , Pm the polygons in the partition produced by Mi U M2. For 
each let ki be the number of flippable diagonals of Pj, and let li be the number of those diagonals 
(among the ki flippable ones) that belong to F. If M2 uses flip'p(M2) flippable edges of T \ Mi, I of 
which are in F, then YT=i = flip(^) " fliPT(^i) - fliPT(^2) and j2T=i k = \F\-l = N/2 -3-1. 
Applying Lemma [52] to each Pi and multiplying the resulting bounds, we obtain a total of at least 
2Eh^T.'^i-T,k triangulations. Hence, 

supp(Mi UM2) > (2/x)S'>x^^'» = (2/^)^/2-3-/^fiip(r)-fiip^{A/i)-fiiPT{M2). 
Next, notice that ^(M2) = 2'-x^^'"PT(^h)~i ^ go ^ave 

2Af/2-3^flip{T)-flipj,{J\/i)-{Ar/2-3) 2^/2-3^KAf-flipT{Jl^i) 

supp(M, U Aft) > = . (13) 



By combining (jl3p with Kasteleyn's method, we obtain 

1 / V- KM2) 



suppfMi U M2) - ^ 2^/2-3x«^-fliPT(Mi) 



N/A 



< 



2N/2-3^KN-f\\pj,{Mi) I 

eeT\Mi 



I E "WM . (14) 



To bound the sum in the parentheses, we write flip2-(Mi) = if + z, where w; = |P fl Mi | and z is the 
number of flippable edges in Mi \ F . Then T \ Mi contains fewer than N /2 — w edges of P, kN — z 
flippable edges not in P, and fewer than 2N — (kN — flip^(Mi)) non-flippable edges. Therefore, 

^ ^ ^lief<^{l^■i2N-{KN-f\\pr{M^)))+x^■{KN-z) + 2^■iN/2-w)). 

eGT\Mi 

Since x < 2, the right-hand side is maximized when w = and z = flipj^(Mi), and it then becomes 
^ J2 Ke? < ^(l'-(2iV-(KiV-flipr(Mi)))+x2.(KiV-flipr(Mi)) + 22.Ar/2) 

eeT\Mi 

= 8 + 2(^2 - 1)k - 2(x^ - 1) • flip7.(Mi)/7V 

2(x2 - 1) flip2.(Mi 



+ 2{x^ - 1)k) 1 



+ 2(x^-1)k N 



j,2 1 flipj(A/i) 
..2 1 \ ,\ „ ITTT^ 



< (8 + 2(2;^ - 1)k) • e N (using l-u< e"" for u > 0) 

= (8 + 2(x2 - 1)k) • (7/x)^-"'P^(*^^)/^ . (15) 
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Combining (jl4p and (jlSp . we get 
1 



< 

J\f]^,i\/2 edge-disjoint 



^ supp(C) ~ ^ supp(MiUM2) 



AfiGA4{T) 

^ MieA^(T) 

To bound the sum in (|16p . we once again use Kasteleyn's technique. This time, we define a weight 
function i^(-) over the edges of T, such that every flippable edge gets a weight of 7, and every other 
edge a weight of 1. Notice that, in this manner, z^(Mi) = j^^'Pri^'h) every Mi € A4{T). We thus 
have 



Mie>i(r) 



\ eer / 



Ar/4 



< (^-(72.f|ip(T) + i.(3Ar-flip(r)))J 

< (^6+-(7'-l)(iV/2 + ^iV)j 

< (6 + 2(72-1)(k + 1/2))^/'. (17) 
Finally, combining (jl6p and ()17p implies the assertion of the lemma. □ 



Note that in the worst case, when k = 0, the bound becomes O ((10 + 27^)^/^). For A; = 0, 
we have 7 = a; • e"*-^^"^-*/"^^, and it is easy to verify that 7 > 1 for 1 < x < 5/4. So the bound is 
actually asymptotically worse than our initial bound of 12^/^, and it continues to be worse when 
K is sufficiently small. As the next subsection shows, in this case the fs-dependent bound from 
Subsection 15.11 becomes small and can be used instead. 



5.3 Integration 

In this subsection we combine the results from the two previous subsections to obtain an improved 
bound for sc(iV). This is done by deriving a connection between 1^3 (T) and flip(r). We start by 
presenting a generalization of Lemma 12.41 

Lemma 5.4 Let c > 1 be a constant such that every set S of an even number of points in the plane 
and a triangular convex hull satisfies sc(5) = O(cl'^l). Then sc{S) = 0{c^^^) also holds for every 
other finite point-set S in the plane. 

Proof. Consider a point set 5. If S has an even number of points, we pick a new point p 
outside the convex hull of S, and denote S' = S U {p}. As mentioned in the proof of Lemma 12.41 
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inserting an additional vertex outside the convex hull of the point set can only increase the number 
of spanning cycles. If S has an odd number of points, we denote S' = S. Notice that, either 
way, S' has an odd number of points. Let Aabc be a large triangle containing S' in its interior, 
and let S" = 5' U {a,b,c}. Again, since inserting an additional vertex outside the convex hull of 
the point set can only increase the number of spanning cycles, we have sc(5') < sc(S" U {a}) < 
sc(S" U {a,b}) < sc{S"). Since S" has an even number of points and a triangular convex hull, 
sc{S) < sc{S') < sc{S") = 0(c^+3) = O(cl^l). a 




Figure 7: Separable edges. 

We also require the notion of separable edges, as presented in [22]. Consider a point set S, 
a triangulation T E T{S), and an interior point p € S. We call an edge e incident to p in T a 
separable edge at p if it can be separated from the other edges incident to p by a line through p. 
An equivalent condition is that the two angles between e and its clockwise and counterclockwise 
neighboring edges (around p) sum up to more than vr. We observe the easy following properties 
(see Figure [7] for an illustration) . 

(50) No edge is separable at both vertices induced by its endpoints. 

(51) If p has degree 3 in T, every edge incident to it is separable (recall that p is an interior point). 

(52) If p has degree at least 4 in T, at most two incident edges can be separable at w. 

(53) If p has degree at least 4 in T and there are two edges separable at p, then they must be 
consecutive in the order around it. 

We are now ready for our main theorem. 

Theorem 5.5 For any set S of N points in the plane, 

sc{S) = O (^10.9247^/^) • tr(5) = O (1.8181^) • tr(5). 

Proof. By Lemma l5.4| we may assume that N is even and that S has a triangular convex hull. 
Recall that 

TeT{S)C€C{T) ' 

We sort the triangulations in the first sum according to the number of interior vertices of degree 3 
that they contain, and get 

(27V+l)/3 

E E E 

i=0 Ter(S) CgC(T) ' 
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(The fact that V3{T) < {2N + l)/3 for every triangulation T is estabhshed, e.g., in [25].) Given a 
triangulation T with V3{T) = i, we can use Lemma ET] to bound '^c&C{T) supp(c) • However, when 
vs{T) is small, the improvement in Lemma [5. II is not significant. In this case we will use instead the 
bound in Lemma 15.31 which, as we now proceed to show, becomes significant when v^iT) is small. 

Consider a triangulation T € T{S). Since S has a triangular convex hull, T contains 3A^ — 9 
interior edges. Notice that an interior edge e is flippable if and only if e is not separable at either 
of its endpoints (this property is equivalent to e being a diagonal of a convex quadrilateral). From 
the above properties of separable edges, we have 

Interior edges Interior vertices of degree 3 Other interior vertices 

flip(r) > -3- V3{T) -2- {N -V3{T) -3) = N -3-V3{T). 



To find for which values of i it is better to use Lemma |5.1|, and for which values it is better to use 
Lemma l5.3( we define t = V3{T)/N and 

_ flip(r) - (iV/2 - 3) ^ (iV-tiV-3)-(jV/2-3) 

^ - N - N " " 

and solve the equation 

/ 3_ / (2-t)(2-t/2) \ "^^^ ^ / (3+(7^-l)(K + l/2))(4 + (x^-l)K) ^ 

I 2* V (1 - i)^ J / V 



where x ~ 1.17965 and 7 = x • e ■4(4+(a:2-i)K) • this will determine the threshold where the two bounds 
coincide. That is, we need to solve the equation 

3 /(2-t)(2-t/2)V-* (3 + (72_i)(^ + i/2))(4 + (x2-1)k) 



2* V (1-0^ 

For this, we use the Wolfram Mathematica software [26], and obtain the solution t ~ 0.1072. 
Moreover, it is easily shown that for i > 0.1072A^ the bound from Lemma l5.ll is smaller, and for 
i < 0.1072A^ the bound from Lemma 15.31 is smaller. In fact, these bounds, in their appropriate 
usage, are all dominated by the common bound for t ~ 0.1072. This, together with (jlSp . imply the 
asserted bound. □ 

By combining Theorem 15.51 with the bound tr(iV) < 30^ [21], we obtain: 
Corollary 5.6 sc(iV) = O (54.5430^). 
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